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Abstract 



The classical Euler-Poinsot case of the rigid body dynamics admits a class of 
simple but non-trivial integrable generalizations, which modify the Poisson equations 
describing the motion of the body in space. These generalizations possess first 
integrals which are polynomial in the angular momenta. 
^^ ' We consider the modified Poisson equations as a system of linear equations with 

^ . elliptic coefficients and show that all the solutions of it are single-valued. By using 

the vector generalization of the Picard theorem, we derive the solutions explicitly 
in terms of sigma functions of the corresponding elliptic curve. The solutions are 
accompanied with a numerical example. We also compare the generalized Poisson 
equations with the classical 3rd order Halphen equation. 

1 Introduction 

As in |lj, we consider the following system 

juj = Jcjj X cjj, j = j X Bu}, a;,7GM'^, (1.1) 
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which is a certain hmit of the Kirchhoff equations describing the motion of a rigid body 
in an ideal fluid. Here u is the angular velocity of the body, 7 is the linear momentum; 
3x3 matrices J, and B are tensors of adjoint masses. The first equation in (|1.1|) are 
just the Euler equations describing the free motion of the body with the inertia tensor 
J. 

In the sequel, J and B are assumed to be arbitrary diagonal matrices. In the special 
case B = Ida, the system becomes the classical integrable Euler- Poinsot case of the rigid 
body motion. In this case vector 7 is a vertical vector fixed in space. Then, in the 
general case, u is solved in terms of elliptic functions, and the 3 independent solutions 
for the vector 7 are elliptic functions and elliptic functions of the second kind, see e.g., 

mm 

Setting M = Jco, the system (jl.ip can be rewritten in the form 

M = M X aM, 7 = 7 X 6M, (1.2) 

where 

a = diag(ai,a2,a3) := J~^, 6 = diag(6i, 62, ^3) := -BJ~^. (1.3) 

It has three independent polynomial first integrals 

Hi:={M,aM), H2:={M,M), Hs:={jn), (1-4) 

Here and below {x,y) denotes the scalar product of vectors x, y G M'^. As the sys- 
tem ()1.2p is divergence free, according to the Euler-Jacobi theorem, for its integrability 
one additional first integral is required. 

In [1] , the authors applied the Kovalevskaya-Painleve method to search for integrable 
cases of the considered system. It was shown that if all the solutions of the system (jl.2p 
are meromorphic, or single- valued, then 

k'^a32ai3a2i + bla32 + blai3 + bla2i = 0, aij = ai-aj, (1.5) 

where k is an odd integer. Geometrically, the above condition describes a quadric in M^ 
with coordinates (61,62,63). In particular, one has b = ka. In this case the equations 
(jl.2p become what can be called the modified Euler-Poinsot system 

M = M X aM, 7 = A; 7 X aM. (1.6) 

As it was shown in [Ij, if the condition (jl.Sp is satisfied, then for odd positive k the 
system (jl.2p possess an additional first integral i?4, which is algebraically independent 
with ()1.4p . It is linear in 7, and of degree k in M, and can be written in the following 
form 

H^ = {P{M),^), (1.7) 

where the vector P{M) is given by 

P{M) = diag(Mi,M2,M3)$fc(M)r, (1.8) 

$fe(M) := (A^^K) . {A^'K) ■ ■ ■ (A'^^i^) . (1.9) 



The matrices K, An are defined as follows 

-na32 &3 -&2 

iC = diag(Mf, M|,M|), An=\ -63 -nai3 61 |, n G N. 

62 -61 -"-021/ 

The constant vector T G M'^ in formula p.Sp spans the kernel of the matrix A].. 
Notice that 

det j4„ = -n(n^a32ai3a2i + ^1032 + b\ai^ + h\a2i). 

As will be shown in Section 8, the case of negative odd k can be reduced to the above 
one. 

In the simplest non-trivial case b = ka with k = 3, we have 

P(M) = P(3) — (pI^\p^^\pPY (1.10) 



with 

pf ) = Ml [(0102 + 8a? - 0302 + a3ai)M^ 

+ (80302 - 80301 + 90i02)M| + (90301 - 30i02 + 30302)M|], 

pP = M2 [(-30302 + 30301 + 9oi02)Mf 

+ (0302 - 0301 + 8O2 + Oi02)M| + (3O3O1 - 30i02 + 90302)M|], 
pP = M3 [(90301 + 30i02 - 80302)Mf 

+ (-3O3O1 + 30i02 + 90302)M| + (03O1 + 8O3 - O1O2 + 0302)M|]. 

In m it was also shown that the vector -P(t) := P{M{t)), where M{t) is a solution of 
the Euler equation in (jl.6p . itself is a meromorphic solution of the Poisson equation in 
()1.6p . Since the solution M{t) in terms of elliptic or hyperbolic functions is well-known, 
the solution P{t) can be found by using (jl.8p . 

In the sequel we will regard the generalized Poisson equations in (II. 6p as a separate 
system of linear equations 
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7 = ^(t)7, A = k\ -a^Ms oiMi e so(3). (1.11 



with the coefficients given by the elliptic functions M{t). (We will not consider the 
special cases when M{t) are hyperbolic functions describing asymptotic motions of the 
Euler top.) 

In the present paper we show that under the condition b = ka {k is odd), all the 
solutions of (II. lip are meromorphic. Our main goal is to give an explicit form of their 
three independent complex solutions in terms of elliptic functions and elliptic functions 
of the second kind (sigma- functions and exponents), as presented in Theorems H] and [6] 
below. 



Additionally, in Theorem [71 we give expressions for the components of the associated 
real orthogonal rotation matrix TZ{t) whose columns satisfy the Poisson equations. 

These equations give rise to a 3rd order ODE for one of the components of the vector 
7. In the final part, we compare this ODE with the best known integrable ODE with 
elliptic coefficients, namely the Halphen equation, and show that, in general, they cannot 
be transformed into eauch other. 

2 General properties of the solutions 

As was already mentioned, the Kovalevskaya-Painleve analysis made in [T] shows that 
for all the solutions of the system (jl.2p or the Poisson equation (jl.lip to be single- 
valued, the condition (|1.5|) must hold and k must be an odd integer. We show that these 
conditions are also sufficient!^. 

Lemma 1. For an arbitrary solution Mi (t) , M2 (t) , M3 (t) of the Euler equations, all 
solutions of generalized Poisson equations (jl.lip are single-valued if and only if k is an 
odd integer and condition ()1.5p is fulfilled. 

Proof. An elliptic or a hyperbolic solution M{t) of the Euler equation has four simple 
poles in the fundamental region. Hence, all singular points of equation (jl.lip on C are 
regular. Since the equation is linear, branching of its solutions can happen only at the 
singular points. 

If k is an odd integer and condition (jl.Sp is satisfied, then all exponents at each 
singular point are integers. However a branching can still occur if the local series solution 
in a neighborhood of a singular point has logarithmic terms. We will show that this never 
happens due to the presence of the first integral (jl.7p . Namely, the integral implies 
that the equation (jl.lip has time dependent first integral 14(^,7) := {P{t),'j), which is 
polynomial of degree k in 7, and P{t) is the corresponding elliptic solution of (jl.lip . 
Assume that P{t) is normalized: {P{t),P{t)) = 1. 

Now take to £ C which does not coincide with a pole of M(t), and a loop s 1— )• t{s) G 
C, s G [0, 1], r(0) = r(l) = to, which encircles once counterclockwise a pole t* of M{t). 
Let r(t) be a fundamental matrix of (jl.lip with the first column proportional to P{t) 
and let T{to) G S0(3, C). Then T{t) G S0(3, C) for all t where it is defined. 

A continuation along the loop r gives a monodromy matrix Mr G S0(3, C): 

T{T{s + l))=T{T{s))Mr. 

For any solution j{t) = T{t)v, v = const G C'^, the integral Ii{t,^) implies 
(P(to),7(to)) = (P(to),r(to)^) = {P{to),T{to)Mr^. 

Since v is arbitrary, this yields P-^(to)r(to) = P'^ {to)^{to)M.T and, due to the ortogo- 
nality of T{t) and the normalization of -P(t), 

(1,0,0) = (1,0,0)X,. 



^In fact, this was already stated in [T], however, without a proof. 



Then, since, Air is also orthogonal, it must have the block structure 

/I 0\ 
Mr=\0 ?? , ^,t9gC, 6'^ + ^'^ = !. (2.1) 
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We now recall that the formal series solution in a neighborhood of the singular point t* 
has logarithmic terms if and only if the monodromy matrix Ai^- is not diagonalizable (see, 
e.g., [2j). However (j2.ip is diagonalizable for any 0,"!? satisfying the above condition. D 

One of the main tools of our subsequent analysis will be a vector extension of the 
known Picard theorem formulated, in particular, in |6l[7]. For our purposes we adopt it 
in the following form. 

Theorem 2. Let Ti and T2 be the common, real and imaginary periods of the elliptic 
solutions Mi(t),Al2{t),M3(t) of the Euler equations. If all the solutions of (II. lip are 
meromorphic, then, apart from the elliptic vector solution ^{t) = P{M{t)) of (jl.lip . 
there exist two elliptic solutions of the second kind ^{t) = G^^'(t), and jit) = G^'^'(t), 
which satisfy 

G'^^\t + Tj) = SjG^^'>{t), G^^\t + Tj) = Sj^G^'^\t), i = l,2, (2.2) 

where 81,82 G C, and, moreover, \8i\ = 1. 

Proof. The existence of at least one vector solution of the second kind, G{t), follows 
from the vector extension of the Picard theorem mentioned above. Let si,S2 be its 
monodromy factors with respect to the periods Ti,T2. 

Let G{t) be another solution of (fTTT]) . and r(t) = (G{t), G{t),P{t)] be a fundamen- 
tal matrix. The monodromy matrices A^i, and A^25 corresponding to periods Ti,T2, 
respectively, are given by 




T{t + Tj) = [sjG, XjG + XjG + p,P, P) = T{t)M„ 
where 



Mj=\ Q y,- 1, .7 = 1,2, 



where Xj^Xj^Pj are certain constants. Observe that, regardless to the values of the 
constants, both monodromy matrices Mi, and M2 are diagonalizable. 

Next, since, by the assumption, all the solutions of (jl.lip are meromorphic, the mon- 
odromy group must be trivial. Therefore, A^i, and M2 commute, and are diagonalizable 
in the same basis. As a result, there exist two independent solutions of the second kind 
GW(i),G(^)(t) forming the fundamental matrix {G'^'^\t),G^'^\t),P{t)). Following the 
general Floquiet theory, the corresponding monodromy matrices A^i,7W2 must satisfy 

det Mj = exp ( / ' Tr A{t)dt] =1, j = 1,2. 



(Here we used the property A{t) S so{3, C).) Hence, since the monodromy of the elhptic 
solution P{t) is trivial, the monodromy factors of G^^^{t), and G^'^'{t) are reciprocal, and 
this imphes ([X^ . 

Further, let for certain constants vi, 1^2 ^ C 

7(t) = zyiG«(t) + 1/2^(2) (i), t£R 

be a real vector solution of the Poisson equation. This means that, for i = 1,2,3, the 
numbers I'lG- (t) and 1^2(^1 (t) are complex conjugated. Then, for the real period Ti, 
the vector ^(t + Ti) is also a real solution. On the other hand, from the above and from 
the monodromy (j2.2p . we deduce that 

7i(t + Ti) = Si iyiGf\t + Ti) + 5f ^ u2Gf\t + Ti), 
which is real if and only if 1 6*1 1 = 1 . D 

3 Algebraic parametrization and elliptic sigma-function so- 
lution for M and P{M). 



We first recall how generic solutions of the Euler equation in (11. 2p can be expressed in 
terms of the Weierstrass sigma functions. We need this fact to derive the general solution 
of the Poisson equation. 

Let us fix a common level of first integrals ()1.4p 

{M,aM) = I, {M,M) = m^, (7,7) = 1 (3.1) 

For a generic values X: ^-^ solutions Mi{t) of the Euler equations are elliptic functions 
related to the elliptic curve E^ given by 

E = {f? = U^{\)], [/4(A):=-(A-ai)(A-a2)(A-a3)(A-c), (3.2) 

where c := Ijw? . Here and below we assume that c 7^ ai, 02, 03. This curve, compactified 
and regularized, has two infinite points oo±. 

A "rational" parametrization of the momenta Mj in terms of the coordinates A, see, 
e.g., [5], have the following form 



Ma = m^- -W— . (3.3) 

y (Oq, - a^)(aQ, - a^) V A-c 

Then, from the Euler equations, we easily deduce that the evolution of A is given by the 
equation 

A = 2mv^-(A - ai)(A - a2)(A - a3)(A - c). (3.4) 

That is, for any A G C, the right hand sides of (j3.3p satisfy the equations ()3.ip . 

For a real motion, i.e., for real values of /, m, and t, if ai < 02 < 03, then one has 
c G (01,03), c / 02. Moreover, 

f[a2,a3], if oi < c < 02, 
I [01,02], if 02 < c < 03. 



The birational map (A,/i) — )■ {z,w), given by 



1 {t2 - 2cTi + 3c^)X + 2cT2 - c^n -3T3 _ fl 

3 A^^^ ' "^ ~ (A - c)2 



transforms the elhptic curve E to its canonical Weierstrass form 

£ = {w^ = Usiz)} , U3{z) := 4{z - ei){z - e2){z - es) = 4z - g2Z - gs (3.6) 

where 



SCa = T2 + CTi — 3(a/3a^ + Caa), ei + 62 + 63 = 0, 

Ti = ai + a2 + as, T2 = aia2 + 0203 + 0301, rs = 010203. 



(3.7) 



The above map sends A = e to 2; = 00, and Oj to 6,, respectively. Then there is the 
following relation between the holomorphic differentials on E and £: 



. dA dz 

1- 



2,/lUX) V4(z-6l)(z-62)(z-63) 

Let us introduce the Abel map 

rp A \ 

^ = W o /7F7TT ' ^'^^'^^ P={\fJ')^E. (3.8) 

Jc zy(74(Aj 

The integrals 

fla := i / — , a = 1,2,3 

A 2yzMA)' 

are the half-periods of the curve E. We choose the sign of the root f/4(A) to ensure 
^1 + 02 + 03 = 0. 

According to (13. Sp . in the case oi < e < 02 < 03 the half-period rJi is imaginary and 
^^2 is real, whereas for oi < 02 < c < 03 the half-period ^^3 is imaginary and ^2 is real. 
In both cases, comparing (j3.8p with (j3.4[) . we get 

u = im{t - to) + n2- (3.9) 

Using the Weierstrass sigma function cr{u) = (7{u\ 2r2i, 2^23), one can write 

A- 6 = const •— — - — — -, where h= / /-— — , (3.10) 

^ = const ■-^--y + ^\ P=r^^. (3.11) 

Moreover, we also have 



= Cq, — -— , a =1,2,3, (3.12) 

A — c a[u) 



where Cq are certain constants and o"a(u) are the sigma- functions obtained from (t{u) 
by shift of u, and by multipHcation by an exponent: 

a„(n):=e'^-" ^^^°~^\ Va = Ci^a), CM = ^^, (3.13) 

where a = 1,2,3. Note that we have 

''^''^^''~m'''~M6'^'^"' ' ^"""^ C7,{0) = a2{0) = asiO) = 1. (3.14) 

see, e.g., |8] or [9]. From ()3.3p . ()3.12p . it fohows that the solutions of the Euler equations 
have the form 

Ma = ha^^^, u = [mt + n2, a = 1,2, 3. (3.15) 

cr(n) 

with certain constants ha which we determine below. In view of (j3.13p . the sigma- 
quotients have the quasiperiodic property 

where 6a,j is the Kronecker symbol. Hence, the coefficients Ma{u) of the Poisson equation 
(jl.lip have common periods 4i7i,4i72. 

Next, using the parametrization p.3p . and the expressions (jl.Sp . for each odd k we 
get the following parametrization for the elliptic solution P{M) = {Pi,P2,Ps) : 



p ^ / {ap - c){a^ - c) X-Ug -F^,»(A) 
^ (aa- ai3){aa- a^)\ A-c(A-c)''' 

where Fg^aW = Po,aY\r=i (-^ ~ /^r.a) is a polynomial of degree s = {k — l)/2, which is 
obtained by substituting (|3.3p into the vector ^kT in (|1.8p and taking the numerator. 
The sum A^ = -Pf (A) + -P|(A) + i-'2^(A) is a constant depending on Ua, and c only. 
In particular, for k = 3, and 6^ = 30^, by using (jl.lOp . we have 

Fii(A) = [3t2 + 4c(c - Ti - 2ai) - 2a2a3]A + c{t2 + 20303) - 4t3 + Sc^oi, 

Fi2(A) = [3t2 + 4c(c - Ti - 202) - 2a3ai]A + c(t2 + 20301) - 4t3 + 80^02, (3.18) 

-F'i3(A) = [3t2 + 4c(c - Ti - 203) - 2aia2]A + c{t2 + 201O2) - 4t3 + 80^03 

and A3 = t| - 4Tir3 + 36cr3 - 48c^T2 + Gic^n. 

Now, applying expressions ()3.1ip . ()3.12p to (j3.17p . we get 

aaju) ^ a{u + Vr,a)cr{u - Vr,a) 
-» o — Cq I I 2/ \ ' [O.iy) 

f^W ;Li ^ («) 

where 

/•Pr,. dA 

frQ = ±i / , , u = i?n,t + S22, 

Jc 2y^lh{X) 

for r = 1, . . . , s. Thus, the components Paiu) have a pole of order k at u = and, like 
Ma{u), they are doubly periodic with common periods 40i,4ri2. We finally have 



Proposition 3. The momentum vector M and the elliptic vector solution P of the 
Poisson equations can be written as 

M -m. ^"(^) P -. ^«(^) rr ^(^ + ^r,a) (t{u - Vr,a) , . 

" - ^^- ^(^) ' ^- - e- ^(^) JUL ^2(^^_^) ^2(^) ' i-^-^uj 

^ (a,/3,7) = (1,2,3), (3.21) 



y^(oQ, - af3){aa - a^) 
where signs of ea are chosen according to the condition 

1 



eie2e3 



-(ai - a2)(a2 - a3)(«3 - ai), 



and n depends on time t via (j3.9p 
Then also 



Paiu) = e„ (p(n) - p(0,)) n (p(u) - p{vi,a,)f , (3.22) 

where p{u) = p{u\g2,g3) is the Weierstrass P-function. 
Here, for any u £ C 

M^{u) + M^{u) + MI{u) = m^, Pl{u) + P|(u) + PI{u) = H, (3.23) 

s 

2 



n = el{p{n^) - p{^p)) J] {p{n^) - p{vi^p)) 

1=1 

s 

+62(p(fi,) - p{n^)) H {p{n^) - p{vi^^)f , (3.24) 

1=1 

for any permutation (a,/3,7) = (1,2,3). 

Remark. According to the rule ()3.16p . the shift u — )• u + 2Qq in the solutions (j3.20p 
is equivalent to flip of signs of some of the constants e^ in such a way that the above 
condition is satisfled. 

Proof of Proposition [21 To calculate the constants ha , Ca in the elliptic solutions (|3.15p , 
(I3.19p . we compare the leading terms of their Laurent expansions near the poles and 
the expansions of the sigma functions. Namely, let to £ C be a pole of the functions 
M(t),P(t), and 5t = t - to. Substituting 

into the equations ()1.6p for M and 7, for any /c G N, one gets 

M(o) G {i(ei, £2, es)^, i(-ei, -£2, £3)^, i(-ei, £2, -£3)^, i(ei, -£2, -£3)^} , (3.25) 
with ea given by ()3.2ip . and P^^' is proportional to M^*^'. 



On the other hand, m view of (j3.14p . near u = we have the expansions 

^ = - + 0(1), 



r=l 



a'^[u) u^ 



with s = {k — l)/2. Since in the above expansions u = \m ■ 6t, comparing them, we 
obtairti 

7 ^a 

Substituting this into (l3T^ . (l3lB . we get dS^Q]). 

The latter, in view of the known relations (see, e.g., [HI [9]) 

c^lju) , . ,f^ . a{u + I3)a{u - (3) 

^^ = Pin) - p{n^), ^2(^)^2(^) = pH - P(/?)' 

implies dS^H). 

Finally, since P{u) is a solution of the Poisson equations, it satisfies the integral 
(j3:23|) . Setting there u = n^ and using ^f322\i one obtains ([SlM]) . D 

Remark. As follows from the formal Laurent solution for M{t) with the coefficients 
(j3.25p . near a pole t = to any vector solution 7(t) = (71,72,73)^ of the Poisson equation 
(jl.lip has the expansion 

const 
7a{t) = j^ {ea + 0{6t)) , a = 1,2,3. (3.26) 

4 Algebraic structure of elliptic solutions of 2nd kind 

Using the algebraic parameterizations ()3.3p and ()3.17p . we obtain 

Theorem 4. 1) If k is a positive odd integer and k > 3, then, apart from the solution 
P{X) in (|3.17p . the Poisson equations (jl.lip has two independent solutions 

^^^^=G^^^ = {G?,G^^\G^i^y, , = 1,2 
which can be represented as the following algebraic functions of the parameter A in (|3.3p . 

daaD 



Ga - ci,„ ^^/^^—-^ exp ^ ^ / W. 



^(,)^ v^OT f_i / 



Co a ^ exp , 



(4.1) 



^ If fact, one can write ha,Ca only in terms of sigma-constants and o"(uj), as it was written for the 
Euler top (the case fc = 1) (see [Hill]), but this process is tedious and requires more calculations. 

10 



Here ci^a,C2,a ore certain constants to he specified below, and 
_ gs+i,a(A)-(A-cr dA 

k 
Qk,aW = ro,a HiX - n,a) = const • (P|(A) + P^{X)) ■ (A - c)'= 

= (a/3 - c){a^ - ap){X - ap)Flp{X) + (a^ - c)(a^ - afi){X - a^)Fl^{X), (4.3) 
(a,/3,7) = (1,2,3), 

where the polynomials Fs^a{X) of degree s = {k — l)/2 are specified in ()3.17p . 

T/ie a6ot;e formula implies that the zeros of the polynomials Qk,a{X) coincide with 
the zeros of Pj (A) + P?;{X). 

2) The differential Wa is a meromorphic differential of the third kind on £ having 
pairs of only simple poles at the points V^^ = {ri^a, iy^tM^V*)") ^ £, i = ^, ■ ■ ■ ,k with 
residues ±1 respectively: 

Resp±W„ = ±l, i = l,...,k. (4.4) 

Finally, qs^i^a{X) in (j4.2p are polynomials of degree s + 1 = {k + l)/2 completely defined 
by the conditions (14. 4p 



The algebraic solutions in the classical case k = 1 will be described separately in 
Section El 

Remark. The polynomials -Fs.ii -^s,2;-^s,3 and Qk,iiQk,2iQk,'i are obtained by the cor- 
responding permutation of ai, a2, 03. Note that their Abel images (u-coordinates) of their 
roots Pr,a,'r'i,a Siie not obtained from each other by the translations by the half-periods 
Qj of the elliptic curve E. 

Proof of Theorem^ 1) According to the kinematic interpretation, the Poisson equations 
in (II. 6p describes the evolution of a fixed in the space vector 7 in a frame rotating with 
the angular velocity ui = BM (also taken in the body frame). Now choose a fixed in 
space ortonormal frame {O, ei, 62, 63}. Let 9, tp, (j) be the Euler nutation, precession, and 
rotation angles associated to this frame so that the corresponding rotation matrix is 

cos (j) cos ip — cos 9 sin ip sin (j) cos (j) sin ip + cos 9 cos ^ sin (j) sin (p sin ( 
TZ = \ — sin (/) cos ip — cos 9 sin tp cos (p — sin (p sin tp + cos 9 cos tp cos (p cos (p sin 9 | . 
sin 9 sin tp — sin 9 cos V' cos 9 

Let -P(t) = (Pi,P2;-P3) be a solution of (II. 6p describing the motion of the vector |P|e3. 
Then, in view of the structure of 7^, and from the Euler kinematic equations, one has 

Pi = |P| sin(/>sin0, P2 = \P\cos(psm9, P3 = |P|cos0, 

and 

■ I |a;iPi+a;2P2 , . _. 

^ = '^' Pf + Pi ' ^'•') 

11 



see e.g., |llj . Hence the thirds components of the other two independent solutions G^^', 
and G^"^' of (jl.lip can be written in the complex form 

G^^'^ = CjSin0exp(±iV'), or G^^'^ = Cj^Pf + P| exp ( ±i f ip dt) , (4.6) 

Cj, and Cj are certain constants, and j = 1, 2. Using the parametrization (j3.17p for Pq,(A), 
as well as relations (13. 3p and (13.41) . we get 



"^ P?(A) + F|(A) 2m.\/CMAy 2 ' ^ ' 

After simplifications this takes the form 

• 1 g.+i,3(A) • (A - cY dX 

lyj at = -- 



2 Qfc,3(A) x/fMA)' 

with 

Q's+i,3(A) = [(a2 - c)(a3 - 02)61 • (A - ai)Fs,i(A) 

m 

+(ai - c)(ai - 03)62 • (A - a2)Fs,2(A)], 
Qfc,3(A) = (^2 - c)(a3 - 02) • (A - ai)Ff i(A) + (ai - c)(ai - 03) • (A - 02)^^ 2(A), 

Pf^ + P| = const ^^S. 

The above implies the formulas (J4.ip ~ (j4.3p for a = 3. Repeating the same geometric 
argumentation for a = 1,2, we get the whole set of formulas of Theorem [H 

2) The differential Wa in (j4.2p has simple poles at A E {?'i,a, • • • ,?'fc,a}, and each of 
them corresponds to two points Vf^^ on E. In view of the degrees of polynomials (5fc,a(A), 
and Qs+i^aj this differential does not have poles at the infinite points cxd-i- on E. Next, 
dX/y^U4{X) does not have neither poles nor zeros on E. Hence Wa is a differential of 
the third kind. 

Next, the residuum conditions (j4.4p are necessary for the solutions (j4.ip to be mero- 
morphic in t, or u and, locally, in A. Namely, let r = A — rj^Q, be a local coordinate on E 
near the root tj^q and the meromorphic differentials have the expansion 

^-=(7 + 0(1)) dr. 

Assume x > 0. Then, as follows from (j4.6p for a = 3, the leading term of the expansion 
of the solution T3 has the form 

const • \/t exp ( — In r j = const • \/t r^' . 

Hence, x must be 1 or 3, 5, ... . Since Tg is an elliptic function of the second kind, 
the total number of its zeros on E must be equal to that of its poles (with multiplicity), 
that is, k, therefore the residuum x must be 1. The same argumentation for a = 1,2 
completes the proof. D 
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5 Sigma- function solutions of 2nd kind 

In order to convert the algebraic solutions of Theorem [¥] to analytic ones, we shall need 
the following formula. 

Proposition 5. Let Kk{\) he a polynomial of odd degree k, and 

•^"-SiiJw- *(A)^.o(A-n)...(A-..) 

he a differential of the third kind on the degree 4 curve E with simple poles at the points 
V, = {rj, ±Y^[/4(rj)), j = 1, . . . ,k with residues ±1 respectively. Let the point (A, /x) G i? 
and u G C be related hy the Ahel map (|3.8p . Then 

r' W = lo^ 'f-^'\- '^; - "^i + 2[C(u;,) + • • • + CM]u + 5u-.i. (5.1) 
7(c,o) cr(^ + wi) • • • a{u + Wk) 

J{\-ri)--- (A-rfc) \/cr(w - wi) • • • a{u - Wk) cr{u + wi) ■ ■ ■ a{u + Wk) 

— , = const -!^ J—— , 

V(A - c)fc tHu) 

(5.2) 
where, as ahove, C{u) is the Weierstrass zeta- function, 6 = K{c)/Q{c), and 

-P- dA 



Wn =\ I — , . (5.3) 

4,0) 2^/ThiX) 



The correct signs of the roots Wj can he chosen from the conditions 
Kk{rj) _ 3i 



{rj - c)Q'^(rj) (ai - c)(a2 - c)(a3 - c) 

The proposition is a reformulation of known relations in the theory of elliptic func- 
tions, its proof is purely technical and given in Appendix 2. 

Note that if the polynomial Kk{X) contains the factor (A — c)^, s > 1, the constant 
6 in (I9.12P is zero. 

Theorem m and Proposition O allow us to formulate the following theorem. 

Theorem 6. 1) The two complex vector elliptic solutions of the second kind of the 
Poisson equations are 

G'"(«)=(g<'>,G<".G<")^, G'^-H") = {Gf\Gf\Gi''y, (5.5) 

GWM^.„e«-n ;'"7''°' , GSHu}=..e-<^-"u t;p"\ , 

fJl o-(^) (^{-^1,0) fJi cr[u) a{wi^a) 



wi,a = ''^ i ;r77fTTf' / = !,..., A;, a = 1,2, 3, 



13 



where ri^a o-f^ the roots of the polynomials QaW in (|4.3p . the signs of Wj^i are defined 
according to (j5.3p . (|5.4p . Next, ta are specified in (j3.2ip . u = imt + ^^2, and 

k k k 

i=i j=i i=i 



Together with ()3.20p . ()3.2ip . i/ie expressions (|5.5p /orm a complete basis of independent 
solutions of the equations (|l.lip . 

2) Let also 

k k k 

^1 = Z] "'J'l' ^2 = X] ^J-2, 5]3 = ^ Wj^z. 

i=i j=i i=i 

The solutions (j5.5p /laiJe i/ie quasi-monodromy 

GW(n + 2f^,) = (-l)^-.,G«(n), Gi2)(n + 2%) = (-l)^-s~i Gi2)(^), 

i = 1,2,3, a = 1,2, 3 
and imply the vector monodromy 

where 



(5.6) 



(5.7) 



si = -exp(2eiJ^i - 2Si?7i) = exp(2e2J^i - 2S2??i) = exp(2e3S7i - 2X3771), 

52 = - exp(2e2fi2 - 2S27?2) = exp(2eif]2 - 2Si??2) = exp(2e3fl2 - 2X3772), (5.8) 

53 = -exp(2e303 - 2X37/3) = exp(2e203 - 2S2773) = exp(2ei03 - 2Si7?3), 

3) If VLj is the imaginary half-period, then \sj\ = 1 . For the real half-period VL2 one has 
\s2\ i^ 1. Moreover, 

J:a-^l3 = ^'y mod{2f]iZ + 2(72Z}, 9^,-6^ = 77^ mod {2771Z + 2772Z}. (5.9) 

/or (a,/3,7) = (1,2,3). 

4) Finally, for any u £ C, 

(G«(7x),G«(tx)) =0, (G(2)(n),G(2)(^)) = 0, 



G«(u)G(,2)(^) = -el Hipiu) - piwr,a)), 



(5.10) 
(5.11) 



and 



r=l 



G«(%)] = sj (-l)i-'-ef n(P(%) - pK«))' 

r=l 
k 

G(,2)(%)]' = sTi(-l)i-^-e?n(p(^-) - p(t.,,«)), 
a,j = 1,2,3. 



(5.12) 
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Moreover, for any n € C, 

k 

-(G«(n),G(2)(^)) = eJHipin,) - p{wi,,)) 

1=1 

k k 

= 4 ll{p{^2) - p{wi,2)) = el llipins) - p{wi,s)) = -U, (5.13) 

1=1 1=1 

where the constant U is defined in (j3.23p . (j3.24p . 

The proof of the Theorem can be found in Appendix 2. 

Remark. One can easily recognize that, for each index a, the components Ga {u), Ga (u) 
have the same structure as solutions of the Lame equation 

— — = (n(n + l)p(u) + 5)A, n G N, 5 = const 
du^ 

with n = k (see, e.g., [12J), namely, 

^1 = n i^rfrrs -p(c(^^)-)) ^ a, = n (4^ eM-c{hs)u)) , 

fj^\a{u)a{hs) J fj^\a{u)a{hs) J 

where the zeros hi, . . . ,h^. satisfy various conditions, in particular, 

p{hi) + --- + p{hk) = kB. 

However, as numerical tests show, the zeros wi^a, ■ ■ ■ -.Wk^a of the solutions (|5.5p do not 
satisfy all the conditions on /ii , . . . , /i^ . Hence Ga {u) , Ga (u) cannot be solutions of the 
Lame equation. 

Thus, if the relation between the Poisson equations (jl.lip and the Lame equation 
(or some of its generalizations) exists, it should be a rather non-trivial one. 

6 The classical case k = 1 

The Poisson equations in this case were first integrated by C. Jacobi [3], who used 
previous results of Legendre (see, e.g., [H]). The case does not fit completely into 
Theorems m [6] because the corresponding meromorphic differentials ()4.2p do not contain 
the factor (A — cY , and the solutions do not have precisely the structure of (15. 5p . 

Namely, now the elliptic solution P is just M{u) given by (I3.20p and the algebraic 
solutions (j4.ip reread 



Gii'2) = Ca^Ml + M^exp U^jWa) 



(a,/3,7) = (1,2,3). 
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Set, for concreteness, a = 3. Using the algebraic parameterization ()3.3p for M(A), from 
(gZD we get 



.aiM2(A) + a2M|(A) 

' M2(A) + M|(A) ' 

(cos — aia2)A + c(aia2 — aia^ — 0203) — 010203 dX 



^^ = '' Mf(A) + M|(l) -^^ 



(6.1) 



(0 + 03 -01 -02)A- (003 -O1O2) 2y^C/4(A)' 

This is a differential of 3rd kind having a pair of simple poles (A*, ^^JU4^{\*)) on E with 

_ ai02 - C03 
+ 03 — 01 — 02 

Observe that in ()6.ip 

(C03 — oi02)A + c(oi02 — 0103 — 0203) — 010203 



(c + 03 — 01 — a2)A — {ca^ — 0102) 
Then, according to Proposition [5l 

'(c,0) 



= 203. 
A=c 



/ '"^ T^3 = log 4^i^^ + 2[C(u;3) + 03]n - vri, 
J(cfi) cr{u + W3) 



Ml + M"^ = const 



^/a{u - W3) a{u + W3) 



cr(n) 

where ^3 is the Abel image of the pole (A*,/i) of W3 with residuum — 1, 

'■^* idA 
W3 



As a result, up to multiplication by —1, the elliptic solutions of 2nd kind are 

^(V), s (j{u — Wa) r... ^ ^ T 

G^^(n) = ea—r-r—, exp (C(^i;a) + aa)u\, 

a{u)a{-Wa) , . 

I ^ (6-2) 

G)^>{u) = ea—r-r—, e^\> [- {Q{w a) + Oajn 

a{u)a{-Wa) 

(compare with (j5.5p ). where Wa denote the Abel image of the pole of the differential Wa 
with the residuum —1, and, as above, u = \mt + il2. 
As follows from item 3 of Theorem [U here 

Wa — Wfi = O^ mod {2^iL + 2^22^}, 
C{wa) - ({wp) =rj-y = Ci^-y) mod {2?7iZ + 2r]2Z}, 
(a,/3,7) = (1,2,3). 

Then, introducing 

w* = Wa — ^a, 0* = ((wa) — Va = ((w*) for any a 

16 



and using the definition (j3.13p of the sigma-functions with indices, one can represent the 
complex solutions (j6.2p in the following form 



G^^H^) = ^^ J."!! /!l exp [{9* + a„)n] 



G^^\r 



a{u) aa{—w' 



(6.3) 



exp[-(6'* + aa)u]. 



By transforming the integrals defining Wa, one can show that 



w 



/(c,o) 2,/lhW 

where oo stands for one of the two infinite points on E. 

Finally notice that being rewritten in terms of theta- functions, the expressions (j6.3 
coincide with the complex solutions presented by C. Jacobi (see also jllj). 



7 Real normalized vector solutions. 



As was shown in Section 2, the elliptic solution P{u) in (j3.20p . (j3.2ip for u = int + Q2, 
t E M, is real. Then, by their construction (see (j4.6p ). the basis elliptic 2nd kind solutions 
G- {u),Gl (u) have opposite arguments: 



Arg(G«(^))=-Arg(Gi2)(t 



(7.1) 



for any a = 1,2,3. Then two independent non-normalized real vector solutions can be 
written in the form 



7(1) (t) = uiG^^\imt + ^2) + 1^2^(2) (imt + O2), 
7(2) (t) = i LiG^^'^iimt + SI2) - U2G^^\imt + ^2) 
with some appropriate constants 1^1 , z^2 • In view of (j7.ip , it is sufficient to set 



iy, = x\G'^^\imt* + n2)r\ 



^2 = X\ 



(2) 



{\mt* + VL2) 



(7.2) 



(7.3) 



for any fixed real t*, any real nonzero Xi ai^d any a G {1, 2, 3}. Then we arrive at 



Theorem 7. A real orthogonal rotation matrix formed by the three independent unit 
vector solutions of the Poisson equations (jl.lip has the form 



7^(i) = -^(7(^)(t),7(')(i),P(t)), 



7«(t) 
7(2) (t) = 



P{t) = Piimt + O2) 
^G(i) (imt + ^2) + yM G^^^ {imt + O2) 



(7.4) 



^^G(i) (imt + O2) + \/^ ^(2) (imt + J^2) 

V-S2 



where P{u) is the elliptic solution (j3.20p . (j3.2ip . and G'-"'^^(n), G^2^(ti) are i/ie elliptic 
solutions of the second kind described in ()5.5p and Theorem iVexi, S2 ^s real and is 
specified in (j5.8p , whereas the constant U is defined in (j3.23p , (j3.24p . 
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Note that the cohimns of TZ{t) form a left- or right-oriented orthonormal basis. 
Proof of Theorem^ Setting in ()7.3p t* = 0,a = 2, and 



S2 



X = y y (P(^2) - P(ai,2)) • • • (p(^2) - p(aA:,2)) , 

in view of (j5.12p , we get z^i = Sg > ^2 = 52 • Then (|7.2j) gives 

7«(0) 
7(2) (0) 
Then, in view of ([530]) and (f5T3|) 



V-S2 



a=l 



a=l 



' = -^(G«(f]2),G(2H^2))=n. (7.5) 



As a result, the real vectors 'y^^' (t) , j^"^' (t) , P (t) all have the same length. By their 
construction, they are all orthogonal. Hence, we obtain the matrix TZ{t) in ()7.4p . D 

8 The case of negative odd k. 

We first note that the case of negative odd k cannot be reduced to the already considered 
case fc > by the trivial substitution t = —T in the Poisson equations in (II. 6p . Indeed, 
this change gives 

^ = -A;7 X aM{-T) 



The elliptic vector solution M{—T) given by (13.20p with u = —irriT + il.2 is neither odd 
nor even, hence one cannot write M{—T) = M{T), and the above equation cannot be 
transformed to the form 

^ = -kjx aM{T). 

Nevertheless, the analysis for positive k is sufficient to cover all the cases. Indeed, 
upon introducing new moments of inertia 



Ai = (Jj + Jk - Ji)Ji, ii,j,k) = (1,2,3), 
the system (jl.6p can be rewritten as 

Auj = —Auj X uj, 'J = kj X uj, 
which, under the change t = —T, gives 

Ao;' = Auj X UJ, 7' = —k'j x uj, (') = -— 

dT 



i.l) 



i.2) 



(8.3) 
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The Euler equations here have the integrals 

(w, Au) = L, {u, A^uj) = M"^ 

with integration constants L,M. Then, applying to these equations the procedure of 
section 3, we express the solutions uj{T) in terms elliptic functions of the curve 

E' = {W^ = -{Z - Ai){Z - A2){Z - A^){Z - C)] 

with the parameters 

A, = 1/A. - "^'"^-"^ 



diCLj -\- fljflfc — CijO-k 



L (r2C-2r3 rg ^ , a2 \ t^ - ATiT-i + At^c 



r2 

3 

Ti = ai + 02 + as, T2 = aia2 + 0203 + 0301, t^ = 010203. 

Here, as above, in (jl.Sp . 

ai = l/Ji, c = l/m^, rri^ = {M,M) = {uj,j'^uj), 

Since u}{T), like M{t), must also be elliptic functions of the original curve E, we get 
the following relation 

d\ 
dt 



2mY^— (X^^a4y(A"^^a^y(A"^^a3y(A"^^c) 
dZ 

~ ~2My^^{X'- Ai){X - A2){X - ^3)(A - C) ' 

Remark. As one may expect, the elliptic curves E,E' with the parameters ai,c and 
Ai, C are birationally equivalent. Indeed, E' is transformed to E by the substitution 

^_ (T2A-2r3)r3 



W = ^ 



rl -4Tir3 + 4T3A' 

(01O2 — O2O3 — Oi03)(ai03 — O2O3 — 0203)(a2a3 — 0103 — O1O2) 



{tI - 4tiT3 + 4 T3 A)2 {tI - 4rir3 + 4 r3 c) 
We stress that the half-periods of E and E' 

/•"" idA , Z'^" idZ 



2yf4(A)' " Jc 2V-(A - ^i)(A - A2){\ - A3)(A - C) ' 
a = 1,2,3 



in general, do not coincide, but only proportional to each other: in view of ([87 
Now comparing (jl.6p or (|8.2p with (|8.3p . we arrive at the following observation. 
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Proposition 8. Let k be an odd negative integer and the vector 

7(T) = j{T\Ai,A2,A3,C,Ai) be a solution of the Poisson equations in (j8.3p with the 
elliptic coefficients oj{T) related to the parameters Ai,C,M. Then 7(i) = 7(— T) is a 
solution of the Poisson equations (jl.lip with the elliptic oefficients Ma{t) related to the 
parameters ai,c,m, and vice versa. 

In other words, the solutions 7(t) of (11. lip with an odd negative k and the parameters 
ai,a2,a3,c, m are given by ^{—t\Ai,A2,A3,C,Ai). The latter are described by the 
formulae of Theorems [4] and [6] corresponding to |A;|, the parameters Ai,A2,A3,C,M, 
and the corresponding roots of the polynomials Fg^aW, Q|fc|,a('^)- 

We stress that although the elliptic vector functions uj{t) for ai, 02, 03, c, m and u;(— t) 
for Ai,A2,A3,C,A4, coincide, this is no more true for the solutions 7(t|ai,a2,a3,c, m) 
and j{-t\Ai,A2,A3,C,M). 

9 Comparison with the Halphen equation 

By using the algebraic parametrization (j3.3p . the generalized Poisson equations (jl.lip 
can be rewritten as 3rd order ODE for one of the components of the vector 7, say 71, 
with the independent variable A G C For general ai,k the explicit expressions for the 
coefficients of the ODE are very long. We give an example for ai = 1, 02 = 2, 03 = 3 and 
k = 3 (c is arbitrary): 

d^ d^ d 

^71 +52 (A) ^71 +51 (A) ^71 + 90 (A)7i =0, (9.1) 

with 

_ 3 3 3 c + 10 1 

^^ ~ 2 (A - 2) "^ 2 (A -3) "^ A-c ~ 10A-15c + 6 + cA ^ A-l' 

3 (138 c -37) 82 c- 119 34 c^ - 276 c + 263 

91 = -TTTT ^T7^ ^ + 



52 (c- 2) (A -2) 24 (c- 3) (A -3) 4 (c - 2) (c - 3) (c - 1) (c - A) 
1 (10 + c)2(10c2 + 2267 c -2666) 



4 (A - 1)2 156 (c - 2) (-3 + c) (7 c - 8) (-15 c + 10 A + 6 + c A) 



5 254 c2- 669 c + 436 

+ 



4(c-A)2 8(7c-8)(c-l)(A-l)' 

27(c-l)(23c-30) 3(c-l)(17c-22) 9 (4c3 - 81 c^ + 171c - 98) 

90 = TTTTT^ :W7 ^TvT + ^TTTTT ^TVT ^vT + 



+ 



104 (A - 2) (c- 2)2 32 (A - 3) (c- 3)2 8 (c - A) (c - 1)2 (c - 2)2 (c - 3)2 
9 (5 c -6) 3(-l + c)(10 + c)5 



16 (c - 1) (-1 + A)2 208 (c - 2)2 (-3 + c)2 (7c - 8) (-15c + 10 A + 6 + c A) 



9 9 (109 c^ - 235 c2 + 106 c + 24) 

^8(c-3)(c-l)(c-A)2 + 32(-l + c)2(7c-8)(A-l) ' 

That is, the coefficients have poles at A = 01,02, 03, c, and at an extra pole defined by 
the condition -15 c + 6 + (10 + c) A = 0. 
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A natural question is how the above 3rd order equation is related with known linear 
equations with elliptic coefficients admitting elliptic solutions of second kind. The best 
known example is the Halphen equation 

^^ + (1 - n2)p(n)|-^ + p'(n)i^* = h^, ^ = ^{u), (9.2) 

where n is integer and h is an arbitrary parameter. As above, p{u) is the Weierstrass 
function. For any such n the 3 independent solutions ^i(n), ^(ti), ^3(u) are elliptic 
functions of 2nd kind with poles of order g = n — 1 at u = 0: 

^„(n) = ^ ^ ^ '\, .^ '-^-^ exp[ (CK 0+- • •+CK"^))n], a = 1, 2, 3. 

The structure of the solutions generalizes that of solutions (|5.5p of our equation ()9.ip . 
So, one can suppose that the equation ()9.ip is a special case of the Halphen equation 
(for /i = 0), when one of its solution is elliptic. 

However, written in the algebraic form with the independent variable z such that 

p(n) = z 



2v/(Z-ei)(Z-e2)(Z-e3) 
the Halphen equation with /i = is 

S £■ d 

A{z - ei){z - e2){z - es)—,^ + g2{z)—^'^ + giiz)—-fi + go{z)^ = 0, 

az"^ dz'^ dz 

5^2 = 362^ + 6ei 62 + 6eie3 + 66263, 5^1 = I22; - 1 - n^, go = -{I - Tn?')/2. 

Hence, the coefficients of the normalized equation have finite poles only at z = ei, 62, 63. 
Taking into account that the equation (jD.ip has 5 poles (which can be reduced to 4 finite 
poles), it cannot be identified with the special case of the Halphen equation. 
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Appendix 1: A numerical example. 

This example was made with Maple by using the functions WeierstrassP (u , g [2] , g [3] ) , 
WeierstrassSigma (u , g [2] , g [3] ) , WeierstrassZeta (u , g [2] , g [3] ) . Consider the sim- 
plest nontrivial case A; = 3 and B = ka. Choose the inertia tensor J = diag(l, 1/2, 1/3), 
i.e., a = diag(l, 2, 3) and l/n? = c = 5/2 € [02, 03], n = 1. Therefore, in the real case the 
parameter AG [ai, 02] = [1; 2]. The elliptic curve E has the form 

/i^ = 2Ui{X) = -2(A - 1)(A - 2)(A - 3)(A - 5/2). 
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The birational transformation ()3.5p . namely, 



1 A + 2 ,_„15^-1 



takes it to the Weierstrass form 

7 10 
w'' = Az^--z + — = A{z + 5/6) (z - 1/6) (z - 2/3), 

with z{\ = 1) = 1/6, z(2) = 2/3, z(3) = -5/6. 

so that the parameters of the Weierstrass functions of E are 52 = 7/3, (73 = —10/27. 
The half-periods aro 

/■"I d\ [^1'^ dz 
Oi = i / — = / — = -1.656638 - 1.415737 • i, 

Jc fJ- Joo W 

/■2/3 J /■- 5/6 J 

02=/ — = 1.656638, JI3 = / — = 1.415737 • i. 

Joo W Joo W 

The corresponding constants r]i = Ci^i) in (|3.13p are 

ryi = -0.4402056 + 0.57199 • i, r/2 = 0.4402056, % = -0.57199 • i. 
This allows to calculate 

aa[u) = ex.p{rjau) — — ^, a = 1,2, 3. 

Next, 

1 _ \/2 _ . _ \/2 

y^(ai -a2)(ai -03) 2 ' ' 2 

From (|3.18|) we have 

Fii{X) = -34A + 167/2, F12 = -48A + 237/2, F13 = -66A + 327/2, 

The Abel images of their zeros on the complex plane u are 

vi = ±0.34497195, V2 = ±0.2898, V3 = ±0.24686 . 

Then the real elliptic solutions for the Euler and the Poisson equations are given by 

a = 1,2,3, u = it + n2, t£R, 

with the indicated above values of the parameters. Note that here |M| = 1 and |P|^ = 
n = 201.062507. 



'in this example most of the float numbers are indicated up to 10 
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Next, the meroniorphic differentials in ()4.2p are 

3\/32l7(2y - 5)(260y2 _ 1295 y + 1613) 



Wi =i 
W2 = i 

^3 ^ 



2(3984 y3 - 29608 y2 + 73363y - 60607) yZ^A) ' 
3\/32T7 (2 y - 5) (560 y^ _ 3136 y + 4331) 

2(23824 y3 _ 194232 y2 + 523569y - 467236) yC/4 (A) ' 
3\/3217(2y-5)(1084y2 -4913y + 5521) 



2(50672 y3 - 356280 y2 + 832461 y - 646139)yC/4(A) 

Tlien, up to constant factors, tlie polynomials Qs,^ in Theorem H] are 

Qg 3(A) = 3167/2 A^ - 44535/4 A^ + 832461/32 A - 646139/32, 
Q3 2(A) = 1489 A^ - 24279/2 A^ + 523569/16 A - 116809/4, 
g3_i(A) = 6723/2 A^ - 99927/4 A^ + 1980801/32 A - 1636389/32 . 

The Abel images of their zeros on the complex li-plane are 

■wi,3 = ±0.21963966, -^2,3 = ±0.309571742, u;3,3 = ±1.232049297, 

wi^2 = H^-i + 0.40885), W2,2.w^^2 = ±(0.272475 ± 0.041941 • i), 

wiA = -0.26888528, u'2,1, ^-3,1 = ±(0.3424574955 ± 0.2549408658 • i). 



Applying the condition (|5.4p . we choose the signs 

u;^ 3 = -0.21963966, '«;2,3 = -0.309571742, 1^3,3 = 1.232049297, 

wi^2 = -0.40885 ± JI3, 'u;2,2 = -0.272475 ± JI3, 

it;3 2 = -0.272475 - J^3, (9.4) 

wi^i = -0.26888528, W2,i = -0.3424574955 - 0.2549408658 • i, 

Wi 1 = -0.3424574955 + 0.2549408658 • i. 



Then 



3 
S3 = ^ Wj^i = 0.7028378946, S2 = -0.95380028 + ^^3, Si = -0.95380028 

and 

3 

03 = ^ C{wj,i) = -7.03775, G2 = -7.477955 - 0.572 • i, 61 = -7.477955 . 
i=i 

Notice that 

S3 - S2 = 1.656638 - 1.41573 • i = Oi = Jli, S2 - Si = Vt^, S3 - Si = O2 (9.5) 

and 

e3-e2 = C(^i), 62 - 61 = C(^3) = r/3, 63 - 61 = c(fi2) = r?2. (9.6) 
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That is, the sums of zeros of (5i(A), (52(A), (53(A) on the complex n-plane differ by the 
half-periods of the curve E, as predicted by item 3 of Theorem [6l Note that for arbitrary 
values of Wij , the relations (|9.5p does not imply (|9.6p . 

The basis complex vector elliptic solutions of the 2nd kind in (j5.5p are 






a{u - wi^a) a{u - W2,a) cr{u - -W3,a) e„« 
" Cr3(u) cr(-U'i,a) cr(-UJ2,a) 0-(-W3,q) 
a{u + Wi^a) a{u + W2,a) CTJU + Wz,a) _0^„ 

a = 1, 2,3, 



with wi a specified in (j9.4p . Next 




{u + 2!^3 



«3 I G« 



-G 



(1) 



S3 = exp(2[eiJl3 - Sir/3]) = -0.962799 + 0.2702178 • i, |s3| = 1, 



'G\ 



(ly 



a 

.G 



(1) 
(1) 



G\ 

{u + 2Vt2) = S2 I -G 



(1) 



(1) 

4'' 



S2 = exp(2[eiJl2 - Si??2]) = 0.402144 • 10' 



-10 



The monodromy of the solutions Gi (u) is inverse to the above one. 



(2), 



Finally, the ortogonal matrix of real solutions is 
1 



n{t) 



(2)/ 



f^^(t)f^Ht)P(t) 



P(i)=P(it + f^2), 



#(t) 



1 



^G«(it + O2) + ^G^^\[t + ^2 

S2 



with ^ = 0.634148 • 10"^ and U = 201.0625 



-^G^^\it + f]2) + \/^G(2)(it + O2) 

V-S2 



Appendix 2: Proofs of Proposition [5] and Theorem [6]. 

Proposition [5] is a reformulation of known relations of the theory of elliptic functions. 
Consider first an elliptic curve £ in the canonical Weierstrass form (|3.6p . 



£ = {w'^ = P^iz) = 4{z - ei){z - e2){z - 63)} , ei + 62 + 63 = 0, 
point P = (z, w) = {z, -y/j^^ij) on it, and the Abel map 

dz 



u 



2y/{z - ei){z - e2){z - 63) 
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: dz. 



(9.7) 



which gives z = p(n). 
Now let 

W 



c(^) 



2^[Wz)Qk{z) 



■ dz. 



Qk{z) = {z- zi)--- {z- Zk), 



(9.8) 



be a meromorphic differential of 3rd kind having pairs of only simple poles at the finite 
points V, 



[Zi 



±2Y^i?3(zj)), i = 1,. . . ,k with residia ±1 respectively. Here qk{2 



bkZ + • • • + &o is a polynomial of degree at most k. 



Theorem 9. If u and P £ £ are related by the map (9.1), then, up to an additive 
constant, 

V = log :S":"^; •••!;" 7"^; + 2[ck) + • • • + c^)]. + --, 0.9) 



a{u + wi) ■ ■ ■a{u + Wk) 



where 



[±Wi) = Zi, - = -p [Wi) 



(9.10) 



C,{u) is the Weierstrass zeta function, p'{u) is the derivative of the Weierstrass P- 
function, and x is the first coefficient in the expansion ofW at the infinite point oo G <S; 
W = {k + 0{u))du, that is, 

qk{z) 



K = lim 

^-^~ Qk[z) 



hh. 



p'{wi), the condition Res^- VF = — 1 is equivalent to 



Proof. In view of 2^JR^i{zi) 
(I91TO . 

It is known ([8]) that the above integral has the form 

^ ^^r 1 aiu — wi)- ■ ■ (t(u — Wk) ^ ^ ^ ^ 

W = \o^—, — ■ — -{ ) — ■ — ^ + Cin + Co, Ci, Co = const. 



a{u + wi) ■ ■ ■a{u + Wk) 

So, it remains to calculate Ci for the differential (19. Sp . Differentiate both parts of (19.9 
by u and evaluate the result at u = (z = oo). Then the right hand side giveqfl 



E 



a'{u — Wi)a{u + Wi) — a'{u + Wi)a{u — w-i) 
a{u — Wi)a{u + Wi) 



+ 2C (wi 



+ x 



u=0 



J2 [C{u - Wi) - C(n + Wi) + 2C{wi)]^^Q + X = X. 



i=l 



Derivation of the left hand side of ([9]) gives 



' d rP=(^''^) \ dz 
dz ,Lr, / du 



~Xz) 



M=0 



lim - , , , 

z-^oo Qk{z) 



*Here we used ("(m) = a'{u)/a{u) and oddness of ^(u) 
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which is precisely bk- □ 

Under a birational transformation {z,w) — )• (A,/i), which sends z = oo to A = c and 



converts £ to the even order curve (I3.2p . 

^^ = [/4(A) = -(A - ai)(A - a2)(A - a3)(A - c), 
the differential (j9.8p takes the form 

Kk{X) dX 



W 



Qfc(A) yC4(A) 



with certain degree k polynomials Kk{\) and QkW = ro(A — ri) • • • (A — r^). Then 
Theorem [9] implies 

^=(^''^) 9,. (A) dX 



(c,o) QfclA) 7CA(A) 

wk= -7=n^dX, 5 = ^, (9.12) 

which is the expression ()5.ip in Proposition El Under the birational transformation ()3.5p . 
the condition (|9.10p takes the form ([5^ . D 

Proof of Theorem [6] 

1). The structm'e of the solutions (j5.5p follows from Theorem S] and Proposition [5j 
Namely, substituting the sigma function expressions (j5.2p . (15.11) for each ^^^^(A) and I¥q, 
into (j4.ip . one obtains these solutions. 

Now notice that, in view of the leading behavior (j3.14p . the Laurent expansions of 
(j5.5p near u = and t = to are 

G(^'^) = ^ + 0(n-+^) = (-l)^'^-^)/^;^^^^ + oat - to)-^-^). 

The leading terms are proportional to those of the required expansions (j3.26p . hence the 

constant factors in the components of Ga ' are correct. 

(1 2) 
The structure of Ga ' {u) implies that they are elliptic functions of the second kind. 

Since all the solutions of Poisson equations are single- valued. Theorem [2] is fully appli- 
cable, hence the vectors G^^''^'{u) must have the vector monodromy (j5.7p with certain 
factors s^. Then 



[u 



gW(u + 2%) = ±5,G«(n), Gi2)(n + 20,) = ±sf G<iHi 
On the other hand, from the quasiperiodicity law of cr{u) we have, in particular, 

G[^\u + 2ni) = siG[^\u), Gf\u + 2ni) = s^^G^^\u), (9.13) 

si = exp(20iOi — 2Si?7i). 
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By the construction of the vectors Ga ' (see (j4.6p ). for any u E C 

g[^\u)Pi{u) + G^^\u)P2iu) + G'i\u)P3iu) = 0. (9.14) 

Here, from ()3.16p and (|3.20p . we have 

P{u + 2(7i) = (Pi(n), -P2{u), -P3{u)f, 
P{u + 2^2) = {-Pi{u),P2iu),-Ps{u)f, 

P{u + 2ns) = {-Pi{u),-P2{u),P3{u)f. 

This, together with the integral (j9.14p . and ()9.13p implies 

G«(n + 2f]i) = {siG^^\u),-siG^^\u),-siGil\u)f. 
Repeating the argumentation for other il.j, we obtain the behavior (j5.6p . ()5.8p . 



3) Let the half-period ilj be imaginary and O2 be real. Then item (2) of Theo- 
rem [2] implies that | s^ | = 1 if we identify the periods Ti , T2 with some of full periods 
4(^1,402,403 of Mi{u). Hence, also \sj\ = 1. By ([SS]), s^ = ±exp(2Gj% - 2X^77^), 
therefore the argument 20jOj — 2Sjr/j is imaginary. Since tjj is imaginary and 772 is real, 
this means that 20^02 — 2T,i'r]2 is real, and S2 is real. 

Next, from ()5.6p . we get, for example, 

-ijy(« + 2O0 = --ijy(n), -^(n + 202) = -ljy(u). (9.15) 

Hence, Gj^ (u)/G3 {u) is an elliptic function with the periods 4^1,202. Its zeros and 
poles in the parallelogram of periods are 

{tfi^i, . . . ,Wk^i,wi^i + 2O1, . . . , Wfc,i + 2ili} and, respectively, 
{w;i,3, . . . , Wk,3,wi,3 + 2O1, . . . , Wk,3 + 2O1}. 

Then, according to the Abel theorem, the difference of their sums must be zero modulo 
the lattice {40iZ + 2^2^}: 

2Si - 2S3 G {40iZ + 2O2Z}, 

therefore. Si - S3 e {20iZ + O2Z}. Further, the case Si - S3 G {20iZ + 2O2Z} is 
not possible (otherwise Gl {u)/G^ (u) would had been a product of an elliptic function 
with the period lattice {2O1Z + 2O2Z} and an exponent, i.e., not doubly-periodic itself). 
Hence, 

Si -S3 = 02 mod {2O1Z + 2O2Z}. 

Applying the same argumentation to the quotients Gg {u)/G^ (""); Gi {u)/G2 (u), we 



arrive at the first part of relations (j5.9p . 
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Next, the quasi-periodicity of a{u) implies 



G\ 



(1) 



G\ 



(1) 



G 



jr{u + 2[7i) = exp [2(Si - Ss)^! - 2(81 - 83)7/1] -f^{u), 



(1) 
3 

.(1) 



G 



(1) 
3 

.(1) 



-^{u + 2^2) = exp [2(Si - E3)l)2 - 2(01 - 63)7/2] -^(^)- 



G 



G 



Comparing this with ()9.15p and using Si — S3 = ^2, as weh as the known Legendre 
relations 



we get the second half of (15. 9p . 

4). Since G^^'{u) is a solution of the Poisson equations, it must satisfy the integral 

G^^^\u) + G^^'^\u) + G'i^\u) = M. 



Then 



G^l^'^iu + 2nf^2) + G^^^'^{u + 2nfi2) + G^^^"^ {u + 2n^2) 



„2n 



.(1) 



(1)^ 



G\^> (u) + G','> (n) + G 



.(1)^ 



M V n G Z. 



Since |s2| 7^ 1, letting above n — )• 00 or n — )• —00, we conclude that M must be zero. 
The same argumentation applied to G^'^>{u), give the second identity in (jS.lOp . 

Let now be a half-period of E, r] = Ci^) ai^d u; G C an arbitrary number. From 
the identity 

a{n + w)a{n-w) 

using the quasiperiodicity of cr{u), we get 



_^2vw_'^^i^-w) 



p(0) - p{w). 



(9.16) 



Applying this to the solutions (|5.5p in the case Q = i^i,w = wn and using (j5.8p gives us 



G^^\n 



^^ 2fT 2C(w„,)n, CT^i^j-W^^i) _ (_-..S^j 2 „.„2S„,y, TT ^^(% "U-^,;) 



k 
1=1 



a'^{Qj)a'^{wa,i) 



fj^a'^{nj)a^{wa,i) 






1=1 



i-lY'^^ etSj ■ {p{nj) - piWaA)) ■ ■ ■ (p(%) - p{Wa,k)), 



.(1), 



.(2), 



i.e., the expressions (|5.12p for G\ {^j)- The expressions for G- (Oj) are obtained in the 
same way. 
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Next, the structure of the solutions (|5.5p gives 
3 3 k 

Y^ G^^\u)G(^\u) = -Y.el n(p(n) - p(^,,„)). (9.17) 

a=l a=l 1=1 

On the other hand, as fohows from (j5.10p and ()5.12p . for any a = 1,2,3 

k k k 

1=1 1=1 1=1 

(a,/3,7) = (1,2,3). 
Then, for it = rij, j = 1, 2, 3, relation ()9.17p yields 

3 k 

- J2 G^aH^j) G^^Hn,) = 2e] n(p(%) - Pi^i,)) ■ 

a=l 1=1 

Since (G''^^(u), G^'^>{u)) is a first integral, we get the relations ()5.13p . 

It remains to prove that the latter equal 11. According to Theorem [5] and relation 

mi, 

k 
Qk,aW = ^o,a HiX - U,a) = const(P|(A) + ^^^(A)) . (A - C)^ 

1=1 

for (a,/3,7) = (1,2,3). This implies, in particular, 

C {p{u) - p{wi,i)) ■ ■ ■ (p(n) - p{wi,k)) = Piiu) + Pi{u) (9.18) 

with a certain constant C. On the other hand, in view of the solutions (I3.22p . one 
obtains. 



Piiu) + Piiu) = el {p{u) - p{^2)) n (P(^) - Pi^i,2)? 

1=1 

s 

+ el (p(n) - p{n,)) H {p{u) - p{vi,3)f . 
1=1 

Letting in the above two relations u — >• oo, taking into account the expansion p{u) = 
l/u"^ + 0(1), and comparing the leading coefficients of l/u"^^, we find 

/o 2 2 2 

^ ^ (02 -ai) (02 -as) (as -ai)(a3 -ai) ^' 

Since Pi{u = Qi) = 0, we have U = P^{u) + P^{u) + P^{u) = if (f^i) + if (f^i). 
Comparing this with (I9.18P for u = $7i, we get 



1=1 
and, therefore, the last equality in (j5.13p . D 
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